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We theoretically study pump-photon-energy-dependent pathways of a photoinduced dimer- 
Mott-insulator-to-metal transition, on the basis of numerical solutions to the time-dependent 
Schrodingcr equation for the exact many-body wave function of a two-dimensional three- 
quarter-filled extended Peierls-Hubbard model. When molecular degrees of freedom inside a 
dimer are utilized, photoexcitation can weaken the effective interaction or increase the den- 
sity of photocarriers. In the organic dimer Mott insulator, K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br, 
the intradimer and the interdimer charge-transfer excitations have broad bands that overlap 
with each other. Even in this disadvantageous situation, the photoinduced conductivity change 
depends largely on the pump photon energy, confirming the two pathways recently observed 
experimentally. The characteristic of each pathway is clarified by calculating the modulation of 
the effective interaction and the number of carriers involved in low-energy optical excitations. 
The pump-photon-energy-dependent pathways are confirmed to be realized from the finding 
that, although the effective interaction is always and slowly weakened, the introduction of 
carriers is sensitive to the pump-photon energy and proceeds much faster. 

KEYWORDS: photoinduced phase transition, metal-insulator transition, dimer Mott insulator, electron- 
phonon interaction 
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1. Introduction 

Nonequilibrium properties of strongly correlated elec- 
tron systems attract much attention. Among them, pho- 
toinduced phase transition dynamics in different groups 
of materials are now deeply and extensively investigated 
from experimental and theoretical aspects. 1 - 1 Their ultra- 
fast and nonequilibrium characteristics are important, 
which we can take advantage of to explore novel func- 
tions. Most of their cooperative characteristics are ex- 
plained on the basis of itinerant electron models. 2 ) 

Photoexcitation energy dependence is always a key 
issue. If the relationship between the density of ab- 
sorbed photons and the amount of photoinduced re- 
flectivity change (such as efficiency and nonlincarity) 
depends largely on the photoexcitation energy, it will 
give a useful hint for optical control in the future. 
For instance, in neutral-ionic transitions in TTF-CA 
(TTF=tctrathiafulvalcnc, CA=chloranil), this relation- 
ship depends on the photoexcitation energy 3, 4 ) and the 
direction of the transition. 5, 6 ) This information con- 
tributes to the knowledge of its transition pathway. 5 ' 7 ) 

If the transition pathway depends on the photoexcita- 
tion energy, it will directly lead to the optical control of 
electronic properties. For instance, Mott insulators are 
known to be generally converted into metals either by 
weakening the effective on-site repulsion or by introduc- 
ing carriers. 8 ^ Photoexcitation is often utilized to intro- 
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duce carriers into one-dimensional Mott insulators and 
induce metallic conductivity. 9, 10 ) One-dimensional Mott 
insulators are, however, special in that infinitesimal on- 
site repulsion leads to the insulating ground state, 11,12 ) 
so that only the introduction of carriers causes metallic 
conductivity. 13 ) 

In two-dimensional Mott insulators, the situation is 
thus drastically different from the above. For instance, 
the Mott insulating ground state in copper oxides is con- 
verted into a superconducting state by chemical dop- 
ing, 8, 14 ) whereas that in organic salts is converted into 
a superconducting state by applying chemical pressure 
to weaken the effective on-site repulsion relative to 
the bandwidth. 15,16 ) Thus, their photoexcitation may 
be able to induce a Mott-insulator-to-metal transition 
via one of these pathways, which generally depends 
on the photoexcitation energy. Indeed, this has been 
suggested to be achieved in deuteratcd «-(d-BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br [BEDT-TTF=bis(ethylenedithio)- 
tetrathiafulvalene] . 17 ) 

In this paper, we discuss how the effective interac- 
tion and the carrier density are modulated by photoex- 
citation in the dimer Mott insulator. It is essential to 
take molecular degrees of freedom inside a dimer into ac- 
count. The organic (BEDT-TTF) 2 X salts basically have 
a three-quarter-filled band. Because of strong dimcriza- 
tion of BEDT-TTF molecules, the K-type salts are gen- 
erally assumed to have the completely filled band orig- 
inating from bonding orbitals and the half-filled band 
originating from antibonding orbitals, which are well sep- 
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Fig. 1. (Color online) Anisotropic triangular lattice for dimer 
Mott insulator. 



arated in energy space. 18 -* Thus, they are often regarded 
as half-filled band systems 19, 20 - ) with a Mott insulating 
ground state. 

It is well known that the effective on-site repulsion 
strength in such a system is given by the transfer inte- 
gral between molecular orbitals inside a dimer in the limit 
of strong on-site repulsion on a molecular orbital. 15 - 1 Its 
modulation by photoexcitation is described in a three- 
quarter- filled-band model, i.e., with two molecular or- 
bitals per dimer. This transfer integral is expected to 
depend sensitively on the distance and the relative ori- 
entation of these molecules. This sensitivity would be 
responsible for the chemical-pressure-temperature phase 
diagram 16 ) in addition to that of interdimer transfer in- 
tegrals. 

The effects of photoexcitation are not so simply de- 
scribed. The resultant charge-transfer (CT) process al- 
ters the stable molecular configuration and thus modi- 
fies the effective on-site repulsion strength. It also gen- 
erally introduces electrons and holes. These two ef- 
fects are always realized, and their relative weights 
depend on the photoexcitation energy. In k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br, both the intradimer and inter- 
dimer CT excitations have broad bands, which overlap to 
a large extent. 21 ' 22 ' Therefore, a single effect is not read- 
ily realized simply by tuning the photoexcitation energy. 

Motivated by this fact, we employ a three-quarter- 
filled extended Pcicrls-Hubbard model on an anisotropic 
triangular lattice. Phonons modulating the effective on- 
site repulsion are treated quantum mechanically, so that 
the intradimer and interdimer CT excitations have over- 
lapped bands even in small clusters. Even under this dis- 
advantageous condition, the selection of the transition 
pathway is shown to be realized in numerical solutions 
to the time-dependent Schrodingcr equation for the exact 
many-electron-phonon wave function. 

2. Three-Quarter-Filled Model for Dimer Mott 
Insulator 

The model we use is the three-quarter-filled extended 
Peierls-Hubbard model on the anisotropic triangular lat- 
tice shown in Fig. 1, 
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where c\ a creates an electron with spin a at site i, 
rii a =(y i(T Ci a , and rii=^2 a . ni a . The parameter U represents 
the on-site Coulomb repulsion, and the intersite Coulomb 
repulsion Vij is assumed to be Vij = Vo/ \ f*j — ry | for 
the four types of pairs (ij) denoted by 61, 62, p, and q 
in Fig. 1. Here, the intermolecular distance is taken from 
ref. 23. The intradimer Vij=Vb\ is the largest among V&i, 
V b 2, V p , and V q . 

The operator b\j creates a phonon, which is as- 
sumed to modulate only the intradimer transfer inte- 
gral tbi- The parameters <?y and u>ij are the corre- 
sponding electron-phonon coupling strength and the bare 
phonon energy, g\,\ and % , respectively. The transfer in- 
tegrals Uj for K-(d-BEDT-TTF) 2 Cu[N(CN) 2 ]Br are esti- 
mated from the extended Huckcl calculation: 24 ) tbi — 
4? - 2t&Mi=-0.2756, t 62 =-0.1047, £ p =-0.1115, and 
i g =0.0404 in units of cV. Here, the transfer integrals 
that are not modulated by phonons are denoted with- 
out the superscript (0). Wc set [7=0.8 and V/,i=0.3. 
Here, we employ the transfer integrals for K-(rf-BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br because this insulating material is 
located in the close vicinity of the metal-insulator phase 
boundary. Because it is impossible to extrapolate from 
the results of the present finite-size system to those in the 
thermodynamic limit, wc cannot judge whether or not 
the present parameter set of tij, U , and Vbi really corre- 
sponds to an insulating ground state in this limit. How- 
ever, the present results and conclusions are not modi- 
fied by the details in these parameters, so that we use 
this parameter set. As for phonons, we take a strong 
electron-phonon coupling gbi=0.06 and a high phonon 
energy w;,i=0.05 to make the intradimer and interdimer 
CT bands overlap to a large extent. 

The time-dependent Schrodingcr equation is solved for 
the many-electron-phonon wave function on the cluster 
of iV=8 sites with a periodic boundary condition (Fig. 1) 
and with the number of phonons restricted to a maxi- 
mum of three at any 61 bond. It is numerically solved by 
expanding the exponential evolution operator with time 
slice dt=0.02 eV _1 to the 15th order and by checking 
the conservation of the norm. This is basically explained 
in ref. 25 for the many-electron wave function, and it is 
now extended to the many-electron-phonon wave func- 
tion | ip(t)). Photoexcitation is introduced through the 
Pcicrls phase, 



c t c . Me/hc)S ir A(t) t c . 
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with Sij = Tj — Ti. The timc-dependent vector potential 
A(t) for a pulse of an oscillating electric field is given by 

F 1 / t 

A(t) = cos(w pmp i)^=^; exp 



pmp 
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where the electric field amplitude F is set parallel to 
the c-axis (the c-component is denoted by F c ), the pulse 
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width is T pmp =35 cV~ 1 =23 fs, and w pmp is the excitation 
energy. 

The transient optical conductivity <r'(u;prb,i) is calcu- 
lated as before, 26,27 ' 



c'(w prb ,i) = - 
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Noj 
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where j = —dH/dA is the current operator, e is a peak- 
broadening parameter set at 0.005, and E = (ip(t) | H \ 

3. Intra- and Interdimer CT Excitations 

Here we discuss the optical modulation of the effective 
on-site repulsion through intradimer and interdimer CT 
excitations. As a first step, we consider an isolated dimer 
consisting of two molecular orbitals overlapped with a 
transfer integral tbi ■ One-hole states consist of the bond- 
ing and antibonding states with energies ± | *6i |. Two- 
hole states consist of three singlet states with energies 
U, (U + V bl )/2 ± v / (^-^i) 2 /4 + 4t 2 1 (~ U, V bl for 
| tbi |"C U, Vbi), and one triplet state with energy Vbi- The 
effective on-site Coulomb energy, defined as E2+Eq — 2Ei 
with E n for the energy of the lowest n-hole state, is given 

by 



U + Vbi 
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4t 2 

^ L bl 



21*61 



(5) 



It becomes [/dim ~ 2 | f&i | +Vbi in the limit of | tbi |<C U, 
Vbi. 

In this strong-coupling limit, the intradimer CT exci- 
tation between the bonding and the antibonding states 
costs 2 | tbi |) while the interdimer CT excitations cost 
2 | tbi | +U and 2 | tbi \ +V b i- Therefore, the interdimer 
CT excitations require higher energies in this limit. How- 
ever, the energy of the intradimer CT excitation is low- 
ered by the second-order perturbation with respect to 
the interdimer transfer integrals, whereas those of the 
interdimer CT excitations are lowered by the first-order 
perturbation. Using realistic values for transfer integrals 
in K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br, we actually find that 
the intradimer CT excitation has a higher energy than 
the interdimer CT excitations, as already assigned in the 
experimental works. 21 ' The optical conductivity is calcu- 
lated by substituting the ground state | ipo) for | ip(t)) in 
eq. (4), as shown in Fig. 2. As clearly shown, the charge 
gap is about 0.18. Because the bare phonon energy is set 
to be Wf,i=0.05, the <r(w) spectrum has peaks at an in- 
terval of about 0.05. The peaks are not exactly located 
at even intervals because they correspond to phonon 
shake-off processes associated with different electronic 
excitations: intradimer and interdimer CT excitations. 
In other words, the intradimer and interdimer CT bands 
are largely overlapped by phonon excitations. In experi- 
mental situations, different modes of phonons with much 
smaller energies contribute to the band broadening, so 
that the spectrum consists of a smooth curve. In the cal- 
culations performed below, the excitation energy uj pl -h is 
set to be an off-resonant energy. As demonstrated be- 
low, the excitations around w=0.3 are mainly due to in- 
terdimer CT excitations, whereas the excitations on the 
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Fig. 2. Optical conductivity with polarization parallel to the c- 
axis. 



high-energy side are mainly due to intradimer CT ex- 
citations. Taking this large overlap between them, the 
picture based on the dimer Mott insulator may not be so 
accurate. At least, the effective on-site repulsion strength 
would deviate from the value in the strong-coupling limit. 

The way in which the effective on-site repulsion is op- 
tically modulated is described rather well in the classical 
picture for phonons. The terms involving phonon opera- 
tors in the Hamiltonian are then approximated to be 



E 

(ij) 
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where sites i and j correspond to two molecules in a 
dimer, the displacement mj is proportional to the classi- 
cal analog of (6^- + &t ), and a^/Kbi corresponds to the 
coupling strength. The force applied to u.y is thus given 

by 

KbiUij 



= a bl yfoLfrg - a \a a l<?) ~ K blU l3 , (7) 
a 

where we rewrite the electron operators as 



bltr,ai a = (Ci 



±c ja 



(8) 



in terms of the bonding and antibonding orbitals in dimer 
I containing molecules i and j. It can be shown, in a 
straightforward manner, that the kinetic term of eq. (1) 
consists of frjj, CT &iCT and a^a^ operators only, whereas the 
current operator j consists of a^bia and b^aia oper- 
ators only. Without Coulomb interactions and without 
photoexcitation, the number of electrons in the bonding 
orbitals and that in the antibonding orbitals were con- 
served. Because the number of electrons in the bonding 
orbitals is larger than that in the antibonding orbitals, 
any photoexcitation reduces So-(^L^ CT — a L a ' CT )> which 
causes force to be applied to Uij in the negative direc- 
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Fig. 3. (Color online) Modulation of effective on-site repulsion 
—At/di m , as a function of the number of absorbed photons 
AE/u pmp , for w pmp =0.30, 0.35, 0.40, 0.45, and 0.50. 



tion, reduces the magnitude of the intradimer transfer, 
and consequently weakens the effective on-site repulsion, 
C^dim- This situation is numerically confirmed below. 

4. Modulation of Effective Interaction and Spec- 
tral Weight after Photoexcitation 

In general, it is difficult to compare numerical results 
with the experimental observation when model calcu- 
lations are based on exact many-electron-phonon wave 
functions, because of the limitation of small sizes of clus- 
ters and small numbers of allowed phonon excitations. 
The present cluster does not show a metal-insulator tran- 
sition because a finite charge gap always exists. However, 
we can evaluate the weakening of the effective interaction 
and the carrier-induced spectral-weight transfer, both of 
which lead to an insulator-to-metal transition in the ther- 
modynamic limit. 

In order to see the modulation of the effective on- 
site repulsion ?7dim in a direct manner, we calculate the 
expectation value of the displacement (6,j + b\j) up to 
£=450. Its maximum decrement, — A(bij + b\-), gives the 
maximum decrement in [/dim, — At/dim- We evaluate it, 



using eq. (5) and the equation tn = t bl — gb\ [bij + b\ 



(0) 
61 

with 4? < 0. We vary the electric field amplitude F c and 
calculate the increment in the total energy AE divided 
by w pmp , which corresponds to the number of absorbed 
photons. We show — AUdim m Fig- 3 as a function of 



AE/u 



pmp 



for different 



y pmp ■ 



As discussed in 8 3, the 



current operator j consists of a\ a bi a and b' ka .ai a oper- 
ators (fc = / for intradimer and fc 7^ I for interdimcr 
CT processes). Both the intradimer and intcrdimer CT 

processes basically reduce X^a-^LA " ~ a L aio ") by two. 
The force applied to phonons is therefore similar between 
these processes. As a consequence, the ratio of — At/dim 
to AE/uj pmp is similar in the range of 0.3 < w pmp < 0.5. 
The effective on-site repulsion is confirmed to be weak- 
ened to a similar extent irrespective of whether charge 
is transferred mainly within a dimer or mainly between 
dimers. 

The number of carriers involved in the optical exci- 
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Fig. 4. (Color online) Increment in conductivity time-avcraged 
and integrated over < Ul < w pr b, A^ajprj,), for (a) cj pr t>=0.01 
and (b) a; pr t, = 0.12, as a function of the number of absorbed 
photons AE/u pmp , for w pmp =0.30, 0.35, 0.40, 0.45, and 0.50. 



tations up to oj pr b is known to be proportional to the 
spectral weight obtained by the integration of the con- 
ductivity over < oj < cj pr b below the charge gap. 
This quantity has been measured, both in equilibrium 14 ) 
and after photoexcitation. 9 ) Here, we calculate the in- 
crement in the conductivity Acr'(u; pr b, i) = <x'(u; pr b,i) — 
c'(w pr b, — 150), time-average it over f50 < t < 750 
Aer'(u; pr b) = (1/600) f 15Q Acr'(u; pr b, t)dt, and integrate 
it over < uj < w pr b, 



iV(w prl 



Aa'(uj)duj 



(9) 



We show iV(a;p r b) in Fig. 4 as a function of AE/u> pmp 
for different w pmp . The energy w pr b is set at 0.01 in 
Fig. 4(a) and at 0.12 in Fig. 4(b) below the charge gap 
of 0.f8 in the ground state. To maintain the numerical 
accuracy, we use weaker electric field amplitudes F c so 
that AE/u> pmp are smaller than those used for the eval- 
uation of — A[/dim- Although the quantity N(uj pi h) in- 
creases with A£/w pmp for any w pmp , its rate depends 
largely on ui pmp . For any u> pr b below the charge gap, 
N(u) pr b) increases rapidly for w pmp =0.3, 0.35, and 0.4 
and very slowly for w pmp =0.5. For w pmp =0.3, 0.35, and 



0.4, the rates arc close to each other. For ui v 



=0.45, the 



rate at w pr b < 0.06 is between those for w pmp =0.3, 0.35, 
and 0.4 and that for a> pmp =0.5 [Fig. 4(a)], whereas the 
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rate at w pr b > 0.07 is close to those for w pmp =0.3, 0.35, 
and 0.4 [Fig. 4(b)]. Namely, the number of carriers in- 
volved in the low-energy optical excitations is increased 
efficiently by 0.3 < w pmp < 0.4, but it is negligibly in- 
creased for aj pmp =0.5. 

This result shows that carriers introduced by photoex- 
citations with cj pmp near 0.3 have low excitation ener- 
gies and are regarded as delocalized. These excitations 
are characterized as interdimer CT excitations. Although 
any CT excitation weakens L7 dim , it requires lattice mo- 
tion and much time. Consequently, if a Mott-insulator- 
to-metal transition is induced, it is mainly through the 
introduction of carriers. On the other hand, carriers in- 
troduced by a photoexcitation with w pmp =0.45 have fi- 
nite excitation energies below the charge gap. A photoex- 
citation with w pmp =0.5 introduces a negligible number of 
carriers. As a consequence, if a Mott-insulator-to-metal 
transition is induced, it is mainly through the weakening 
of Udim- This excitation is characterized as an intradimcr 
CT excitation. The characteristic of such a CT excitation 
seems to vary continuously as a function of the photoex- 
citation energy. 

5. Conclusions 

The present theoretical study is motivated by the 
photoinduced Mott-insulator-to-metal transition ob- 
served in the organic dimer Mott insulator, k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br. 17 ) Using a two-dimensional three- 
quarter-filled extended Peierls-Hubbard model, we con- 
sider excitation-energy-dependent transition pathways. 
Bearing the above material in mind, we introduce quan- 
tum phonons that modulate intradimcr transfer integrals 
so that the intradimcr and the interdimer charge-transfer 
excitations have broad bands that overlap to a large 
extent. The transient quantities are obtained from the 
numerical solution to the time-dependent Schrodingcr 
equation for the exact many-electron-phonon wave func- 
tion on a small cluster. They indeed depend strongly on 
the excitation energy. 

The spectral-weight analysis shows that the number 
of delocalized carriers is increased efficiently when the 
excitation energy is around the peak in the conductivity 
spectrum. If the system were large enough, this would 
lead to a photoinduced Mott-insulator-to-metal transi- 
tion mainly through band-filling control. On the other 
hand, when the excitation energy is away from the peak 
on the high-energy side, few delocalized carriers are intro- 
duced. Slow lattice motion is induced by any CT excita- 
tion and modifies intradimer transfer integrals and weak- 
ens the effective interaction. Consequently, if a Mott- 
insulator-to-metal transition is induced by such a pho- 
toexcitation, it is mainly through bandwidth (relative 
to the effective interaction strength) control. Namely, 
photoexcitation-energy-dependent pathways are realized 
from the fact that, although the effective interaction is 
always and slowly weakened, the introduction of carriers 
is sensitive to the photoexcitation energy and proceeds 



much faster. 
Acknowledgment 

This work was supported by Grants-in-Aid for Sci- 
entific Research (C) (Grant No. 19540381 and No. 
23540426), Scientific Research (B) (Grant No. 20340101) 
and Scientific Research (A) (Grant No. 23244062), 
and by "Grand Challenges in Next-Generation Inte- 
grated Nanoscience" from the Ministry of Education, 
Culture, Sports, Science and Technology of Japan, 
and the NINS program for cross-disciplinary study 
(NIFS10KEIN0160). 



18 
19 

20 

21 

22 
23 
21 
25 
26 
27 



Y. Tokura: J. Phys. Soc. Jpn. 75 (2006) 011001. 

K. Yoncmitsu and K. Nasu: J. Phys. Soc. Jpn. 75 (2006) 

011008. 

T. Suzuki, T. Sakamaki, K. Tanimura, S. Koshihara, and 
Y. Tokura: Phys. Rev. B 60 (1999) 6191. 

P. Huai, H. Zheng, and K. Nasu: J. Phys. Soc. Jpn. 69 (2000) 
1788. 

H. Okamoto, Y. Ishigc, S. Tanaka, H. Kishida, S. Iwai, and 

Y. Tokura: Phys. Rev. B 70 (2004) 165202. 

K. Yoncmitsu: Phys. Rev. B 73 (2006) 155120. 

H. Uemura and H. Okamoto: Phys. Rev. Lett. 105 (2010) 

258302. 

M. Imada, A. Fujimori, and Y. Tokura: Rev. Mod. Phys. 70 

(1998) 1039. 

S. Iwai, M. Ono, A. Macda, H. Matsuzaki, H. Kishida, 
H. Okamoto, and Y. Tokura: Phys. Rev. Lett. 91 (2003) 
057401. 

H. Okamoto, H. Matsuzaki, T. Wakabayashi, Y. Takahashi, 

and T. Hasegawa: Phys. Rev. Lett. 98 (2007) 037401. 

J. Solyom: Adv. Phys. 28 (1979) 201. 

H. J. Schulz: Int. J. Mod. Phys. B 5 (1991) 57. 

N. Maeshima and K. Yonemitsu: J. Phys. Soc. Jpn. 74 (2005) 

2671. 

S. Uchida, T. Ido, H. Takagi, T. Arima, Y. Tokura, and 

S. Tajima: Phys. Rev. B 43 (1991) 7942. 

K. Kanoda: Hyperfine Interact. 104 (1997) 235. 

K. Kanoda: J. Phys. Soc. Jpn. 75 (2006) 051007. 

Y. Kawakami, S. Iwai, T. Fukatsu, M. Miura, N. Yoneyama, 

T. Sasaki, and N. Kobayashi: Phys. Rev. Lett. 103 (2009) 

066403. 

H. Kino and H. Fukuyama: J. Phys. Soc. Jpn. 65 (1996) 2158. 

H. C. Kandpal, I. Opahle, Y.-Z. Zhang, H. O. Jeschke, and 

R. Valcntf: Phys. Rev. Lett. 103 (2009) 067004. 

K. Nakamura, Y. Yoshimoto, T. Kosugi, R. Arita, and 

M. Imada: J. Phys. Soc. Jpn. 78 (2009) 083710. 

D. Faltermeier, J. Barz, M. Dumm, M. Dressel, N. Drichko, 

B. Petrov, V. Semkin, R. Vlasova, C. Meziere, and P. Batail: 

Phys. Rev. B 76 (2007) 165113. 

M. Dumm, D. Faltermeier, N. Drichko, M. Dressel, C. Meziere, 

and P. Batail: Phys. Rev. B 79 (2009) 195106. 

T. Mori, H. Mori, and S. Tanaka: Bull. Chcm. Soc. Jpn. 72 

(1999) 179. 

M. Watanabe, Y. Nogami, K. Oshima, H. Ito, T. Ishiguro, 

and G. Saito: Synth. Met. 103 (1999) 1909. 

K. Yonemitsu and N. Maeshima: Phys. Rev. B 79 (2009) 

125118. 

H. Matsucda and S. Ishihara: J. Phys. Soc. Jpn. 76 (2007) 
083703. 

K. Onda, S. Ogihara, K. Yonemitsu, N. Maeshima, 
T. Ishikawa, Y. Okimoto, X.-F. Shao, Y. Nakano, H. Yamochi, 
G. Saito, and S. Koshihara: Phys. Rev. Lett. 101 (2008) 
067403. 



